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^ ■ Abstract 

o, 

O , Closed-form expressions are obtained for the generating function of 

close-packed dimers on a 2M x 2N simple quartic lattice embedded 
pq I on a Mobius strip and a Klein bottle. Finite-size corrections are also 

^ ■ analyzed and compared with those under cylindrical and free boundary 

'nJ" I conditions. Particularly, it is found that, for large lattices of the same 

size and with a square symmetry, the number of dimer configurations 
on a Mobius strip is 70.2% of that on a cylinder. We also establish 
^^ ■ two identities relating dimer generating functions for Mobius strips and 

Q^ ' cylinders. 
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1 Introduction 

It is well-known that the problem of closed-packed dimers on planar graphs 
can be solved, meaning that its generating function can be formulated as a 
Pfaffian, the square root of an antisymmetric determinant ||l|. This leads 
to the possibility that the generating function can be evaluated in a closed 
form for lattices with periodic structures. Indeed, Kasteleyn [|l| has suc- 
ceeded in obtaining the closed-form expressions of the generating function 
for the simple quartic lattice with both free and toroidal boundary con- 
ditions, and McCoy and T. T. Wu [^ extended the result to cylindrical 
boundary conditions. Temperley and Fisher Q, and Fisher [Q] have also 
solved the free boundary case independently, although the generalization of 
the field-theoretic approach used in Q to general graphs is less apparent. 

While the solutions under various boundary conditions yield the same 
per-site generating function in the bulk limit, their finite-size corrections 
are nevertheless different. In this context there exists two unique surface 
structures, namely, the Mobius strip and the Klein bottle, which have es- 
caped heretofore attention. The Mobius strip is a non-orientable surface 
which has a single side and a single boundary edge, and the Klein bottle, 
also non-orientable, has a single side and no edges. In view of the increasing 
interest on finite-size corrections, particularly their interplay with the con- 
formal field theory H, it is of pertinent interest to study lattice models on 
these surfaces. In this Letter we consider this problem. We solve exactly the 
dimer statistics for a simple quartic net embedded on a Mobius strip and 
on a Klein bottle. We also analyze its finite-size corrections and establish 
two identities relating generating functions under the Mobius and cylindrical 
boundary conditions. 

2 Formulation 

Consider first a 2M x 2N simple quartic lattice C embedded on a Mobius 
strip. Here, for convenience, we assume the number of sites be even in each 
direction. The Mobius strip of length 2N and width 2M is shown in Fig. 1, 
where repeated sites are denoted by open circles. We consider close-packed 
dimers on C. 

Let the dimer weights be zi and Z2, respectively, for dimers in the hor- 
izontal and vertical directions. We need to evaluate the dimer generating 
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Figure 1: A Mobius strip of a simple quartic net of length 2N and width 2M. 
Sites denoted by open circles are repeated as indicated by the labeling. The 
edge orientation satisfies the Kasteleyn clockwise-odd rule for superposition 
polygons, and rectangles of dotted lines denote dimer cities. 
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where n^ and n„ are, respectively, the numbers of horizontal and vertical 
dimers subject to the sum rule n/j + n^, = 2MN . The summation in ([l|) is 
taken over all close-packed dimer configurations. 

A key element in the Pfaffian formulation of the generating function 
is to orient edges of C such that every superposition polygon, formed by 
superimposing two dimer configurations, must contain an odd number of 
arrows in the clockwise direction, the clockwise-odd rule. Kasteleyn Q has 
shown that when this rule is satisfied, then the generating function of close- 
packed dimers can be written as a Pfaffian. For the Mobius strip £, it 
can be shown that the clockwise-odd rule is satisfied if we adopt the edge 
orientations shown in Fig. 1. Then, we have the result 
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izi,Z2) = PiiA) = Jdet\A\ 
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where A is a AMN x AMN antisymmetric matrix whose elements can be 
read off from Fig. 1. 

To see the general structure of A, it is instructive to first write down the 
matrix A in a simpler case. For this purpose consider the 2x4 Mobius strip 



(M = 1, A^ = 2) shown in Fig. 2. Reading off from Fig. 2, we obtain 
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a result which can be verified by explicit enumerations. 
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Figure 2: A 2 x 4 Mobius strip. 
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To write down A for general M and A^, we adopt "dimer cities" formed 
by pairing sites {2n, 2n + 1}, n = 1, 2, • • • , 2MN as indicated in Fig. 1, and 
introduce 2x2 matrices 
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Further introduce N x N matrices 
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and the 2M x 2M matrix 
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J2M = 
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Then, the matrix A can be written in terms of a hnear combination of direct 
products of smaller ones, 



A = A2N O hM + a{0, 1) ® /at® 



F2M 



^2M 



+ B2N®J2M, (8) 



where 



B2N 



a(0, 0) /tv + a(l, 0) F^v + a(-l, 0) F^ 
-a(l, 0) i^TV + a(-l, 0) ® K%, 



(9) 



F^ and KJj are transposes of F/v and Kat, and /at is the N x N identity 
matrix. Note that the first two terms in (^) are the corresponding A matrix 
for free boundary conditions. 

We next consider a Klein bottle constructed by connecting the upper 
and lower edges of the Mobius strip of Fig. 1 in a periodic fashion with 2N 
extra vertical edges. The clockwise-odd rule again holds if these edges are 
oriented inversely to those in the same column. Thus, we have the matrix 
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A - a(0, 1)<E)In<S) {K2M - K^m)- 



(10) 



3 Evaluation of determinants 

We first evaluate det|A|. Since J2M commutes with F2M — F^M^ ^^'^ ^^ ^ 
result they can be diagonalized by applying a common similarity transfor- 
mation. Introducing the 2M x 2M matrix U whose elements are 
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we find 

{U l2MU)m,m' = Sm,m' 
(U~ J2MU)m,m' = i ( — 1) ™ 5m,m' 
{U{F2M-F^M)U~^)m,m' = {2icOS ^m) &m,m' , m = 1, 2, ..., 2M,(12) 

where 

rnvr 

"^^ = (2mTT) • ^''^ 

Thus, we can replace the 2M x 2M matrices in (|8|) by their respective 
eigenvalues, and express det|^| as a product of the replaced determinants, 
namely. 
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Introducing (0) and the N x N matrix 



(14) 



'^N 



FN+i{-iy 





/ 







1 


•• 


■ °\ 














1 •• 


• 


M+m+lj^^ = 


























•• 


• 1 




\ii- 


_-^\M+m+l 





•• 


• 0/ 



(15) 



we can rewrite (|Tj) as 
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det\A\ = Yl det 
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[a(0, 0)+(2i cos (I)m)ai0, l)]«)I^+a(l, 0)®r^+o(-l, 0)8)T^ 

(16) 
where T^ is the Hermitian conjugate of T/v- 

Now Tjv and T^ commute and they can be simultaneously diagonalized 
with respective eigenvalues e " and e~ ", where 
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It follows that we have 



M+m+l 



(4n - l)7r 

2iV 



n = l,2,...,iV. 



(17) 



zMob 



(2M,27V)Ul'^2j 



detl^l 

2M N 

nn 

m=l n=l 



det 



a(O,O) + (2icos0m)a(0, 1) 



+a(l,0)e*^" +a(-l,0)e- 



M N r 

nn 

m=l n=l 



Azi sin 



2 (4 n - l)7r 

4A^ 



+ 4^2 cos 



mvr 



2M + 1 



1/2 



(18) 



where we have used cos^ (pm = cos^ (j)2M-m+i to effect the square root. It is 
readily verified that ( |l8|) reduces to (Q) when M = 1, N = 2. 

The determinant detl^^'*^*"! of the matrix ( p!o[ ) can be evaluated in a 
similar fashion from which we obtain 
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As a comparison, the generating functions for 1M x 2N lattices with free 
boundaries O and cylindrical boundary conditions (periodic in the N direc- 
tion) [H are, respectively. 
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The similarity between expressions (|l8| ) - (20) is striking, with the only 
difference being the trigonometric factors which can be associated to the 
respective boundary conditions. We remark that for (2M — 1) x 2N lattices 
we have 101 
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These expressions give rise to the same bulk "free energy" 
/bulk = , lip (4MiV)"^lnZ(2M,2Af)(^i,^;2) 

= -^ I ^ d9 [ ^ d(t) ln(4z^ sin^ + 4z^ cos^ 
IQii'^ JO Jo 



(21) 



(22) 



Particularly, the number of ways Wgxs that an 8 x 8 checkerboard can be 
covered by 32 dominoes is obtained by setting M = N = A and zi = 2:2 = 1; 
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yielding 

VF^^^Ij' = 12 988 816 = 2^^ X (901)^ 

Wl^°^ = 46 029 729 = 47 x 271 x 3617 

P^g^^ = 71 385 601 = (8449)^ 

Wl^^^^"^ = 220 581 904 = 2^ x (3713)^ 

^t^oidai ^ 3^^ 353 3^2 = 2* X 3^ X 135353, (23) 

where the number for toroidal boundary conditions is computed using the 
result of [|. 

Two identities: There exist curious identities connecting generating func- 



tions for Mobius strips and cylinders. Brankov and Priezzhev [10| have 
proposed the identity 



-^(2M,4Ar)(^l'^2) 



yMob I \ 



(24) 



which they established up to the order of C3 in the large M, N expansion (|26 
The identity (^4|) can now be rigorously established using (|l^) and (pO|) . The 
proof follows from the fact that in the product over n = 1, 2, • • • , \N on the 
LHS, the factor sin^[(2n - l)7r/4Af] reproduces the set sin^[(4n' - l)7r/4iV], 
n' = 1, 2, • • • , A^, exactly twice. For (2M — 1) x 2N lattices we have instead 
the identity (see Note added below) 



^(2Af-l,4Ar)(^l'^2) — „ 



'^(2M-l,2Af)(^l'^2) 



(25) 



4 Finite-size corrections 

For large M and N, we expect to have 

In ^(2M,2W) = 4MiV/buik + 2Nci + 2Mc2 + C3 + 0{1/N), (26) 

where ci, C2 and C3 are constants dependent at most on 

r = Z2/Z1 and ^ = N/M, (27) 

the aspect ratio. The computation of finite-size corrections for products of 
the form of ( \lq j is standard Q, |^, 0], and is outlined in the following. 



We first take out a factor zf in ( [l8| ) and carry out the product over 
n by introducing the identity H 
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with a = ■k/2N . Tliis leads to 
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It follows that we have 
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(31) 
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1 + x-^^^(t, 



(32) 

We next carry out the first summation in ( |3^ using the Euler-MacLaurin 
summation formula 
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with a = and 6 = 7r/(2M+ 1). The first term and the leading contribution 
of the second term lead to the following expressions for the bulk free energy. 
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where the second line is written down by symmetry. 
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After some manipulation following |^] , this leads to the following for the 
Mobius boundary condition: 
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where ??cj is the theta function 



^3(u, g) = 1 + 2 ^ g"' cos(2nu), 



(36) 



n=l 



with q = e"'^^'^' ^. Note that C2 = is a consequence of the fact that there 
is no boundary edge in the M direction. 

For the cylindrical boundary condition periodic in the N direction, one 
has in place of (|29|) , 
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from which one obtains 
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It follows that, for large M, N, we have 
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Particularly, for M = N and zi = Z2 = 1, we have q = e"'^'^ and (^) gives 
the ratio of the numbers of dimer configurations under the two respective 
boundary conditions as 0.701845. 

For free boundary conditions we obtain 
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Here, use has been made of the second line of p4|). These expressions 
reproduce those obtained by Fisher [Q] and Ferdinand Q . Similarly we have 
for the Klein bottle 

cf^*"(r,e) = cf°^(r,0+vrTC/8. (41) 

If one takes the limit of A^ ^ cxd (M -^ oo) first while keeping M (N) 
finite [^, one obtains from (26) 



lim — -Z(2M,2^) = 2Mfbuik + ci + Ai/2M + 0{l/M^), 
hm -i-z = 2Nfhuik + C2 + A2/2N + 0{l/N^). (42) 

Af — >oo ZlVI 

One finds after some algebra the results 
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where we have also included results for toroidal boundary conditions worked 
out from results of M. The tabulated values of Ai and A2 are consistent 
with the finding of [^ that the periodical boundary condition (in the lateral 
direction of an infinite strip) is distinct from the free or fixed, now including 
the Mobius, boundary conditions. While it is tempting to extract the central 
charge c from Ai and A2 resulting in c = r and r~^, respectively, we remark 
that the extraction is superfluous since results of Q apply to systems at 
criticality while the present dimer system does not exhibit a critical point. 
Note added: After the submission of this paper we learned of a recent 
preprint by Tesler 1^] in which the dimer generating function Z^J^^j^Jzi, Z2) 
is enumerated for all M,M (even or odd) as linear combinations of two 
Pfaffians with the results given in terms of a g-analogue of the Fibonacci 
number. We have verified that for M,M even Tesler's result is the same as 
our expression (^). For Ai odd and M even, Tesler found the solution also 
given by (|l|) but with an extra factor z^^^"^ with 2M ^ M, 2N ^ M, 



and the product over m taken from 1 to (A4 + l)/2. We remark that this 
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leads to the identity ( p5[ ) upon using (^). For Ai even and M odd, however, 
the solution assumes a more complicated expression Q. The Mobius strip 
problem has also been studied by Brankov and Priezzhev pH] in the context 
of the large M, N expansion (|2q) of the free energy. 
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